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Abstract 



The production of or Z bosons in association with two jets is an important 
background to the Higgs boson search in vector-boson fusion at the LHC. The purely 
electroweak component of this background is dominated by vector-boson fusion, which 
exhibits kinematic distributions very similar to the Higgs boson signal. We consider 
the next-to-leading order QCD corrections to the electroweak production of ivijj and 
events at the LHC, within typical vector-boson fusion cuts. We show that 
the QCD corrections are modest, increasing the total cross sections by about 10%. 
Remaining scale uncertainties are below 2%. A fully- flexible next-to- leading order 
partonic Monte Carlo program allows to demonstrate these features for cross sections 
within typical vector-boson-fusion acceptance cuts. Modest corrections are also found 
for distributions. 



1 Introduction 



Vector-boson fusion (VBF) processes have emerged as a particularly interesting class of 
scattering events from which one hopes to gain insight into the dynamics of electroweak 
symmetry breaking. The most prominent example is Higgs boson production via VBF, that 
is, the process qq qqH, which can be viewed as quark scattering via t-channel exchange of 
a weak boson, with the Higgs boson radiated off the W or Z propagator. Alternatively, one 
may view this process as two weak bosons fusing to form the Higgs boson. The kinematic 
characteristics of this process are very distinctive: two jets, in the forward and backward 
region of rapidity, with the Higgs boson decay products in the central region of the detector. 
This characteristic signature greatly helps to distinguish these Hjj events from backgrounds. 
Higgs boson production via VBF has been studied intensively as a tool for Higgs boson 
discovery P |2] and the measurement of Higgs boson couplings |^ in pp collisions at the 
CERN Large Hadron Collider (LHC). 

Analogous to Higgs boson production via VBF, the electroweak production of a or 
Z plus two jets, with the requirement that the weak boson is centrally produced and that 
the two jets are well separated in rapidity, will proceed with sizable cross section at the 
LHC"^. The decay leptons inW —>■ iu^ and Z — > lead to the final states iiy^jj and 
{i = e, fi, t). These processes have already been considered in the literature at leading order 
(LO). To name but a few examples, they have been studied in the investigation of rapidity 
gaps at hadron colliders [3 13 El , as a probe of anomalous triple-gauge-boson couplings 
or as a background to Higgs boson searches in VBF fUl IQ- In this last case, the 
d-Vnii final state with an unidentified charged lepton, or vivijj events from Z^viVi decay, 
form a background to invisible Higgs boson decay (see e.g. Ref. j^). T'^T'jj events are a 
background to the decay H —>■ t^t~ ^U], and also to H W^W~ when the W^s and the r's 
decay leptonically ^1]. In these examples, off-shell corrections to Z ^t^t^ decay need to 
be included, since a Higgs boson mass in the range 114 GeV < mu < 200 GeV, well above 
the Z peak, is favored by electroweak data [TB] . 

While a LO analysis is perfectly adequate for exploratory investigation, precision mea- 
surements at the LHC require comparison with cross-section predictions which include 
higher-order QCD corrections. A poignant example is the extraction of Higgs boson cou- 
plings, where expected accuracies of the order of 10%, or even better |3|, clearly require 
knowledge of the next-to-leading order (NLO) QCD corrections. In addition, one would like 
to exploit W and Z production, in VBF configurations, as calibration processes for Higgs 
boson production via VBF, namely as a tool to understand the tagging of forward jets or 

^Another source of Wjj or Zjj events are QCD processes at order a^a, soraetimes called QCD Vjj 
production. Within typical VBF cuts, cross sections for these QCD processes are only somewhat larger than 
those for electroweak production 4 . One thus needs to calculate NLO QCD corrections for both sources 
independently, and as a function of phase space. For the QCD processes this was done in Ref. 
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the distribution and veto of additional central jets in VBF (see e.g. Ref. [HE])- In fact, these 
processes share the same color structure: two colored quarks are scattered via the exchange 
of a colorless boson in the t-channel. The pattern of soft gluon radiation is then the same. 
Understanding the gap-survival probability in the known case of W and Z production can 
give insight on the gap survival for the case of Higgs boson production. The precision needed 
for Higgs boson studies and for the understanding of radiation patterns then requires the 
knowledge of NLO QCD corrections for Wjj and Zjj production as well. 

The NLO QCD corrections to the total Hjj cross section from VBF has been known for 
many years jHj. In a recent paper we presented the calculation of these corrections in 
the form of a fuUy-fiexible parton-level Monte Carlo program which allows the determination 
of NLO corrections to arbitrary (infrared-safe) distributions. Here, we extend this work and 
describe the calculation and first results for such corrections to Wjj and Zjj production in 
VBF configurations. To be precise, since the decaying weak bosons are spin-one particles, 
in order to retain all the possible angular correlations between the final state particles, 
we consider the electroweak processes pp i^u^jjX and pp^i^i~jjX at NLO. At LO, 
Feynman graphs for one such process, uc^dcW'^ ,W'^ l^vi, are shown in Fig. Using 
the terminology introduced in we consider bremsstrahlung (a, b, c), fusion (d) and 
multiperipheral (e, f) diagrams. We neglect diagrams corresponding to conversion, abelian 
and non-abelian annihilation, since these qq annihilation contributions are negligible when 
we impose VBF cuts, as explained in detail in Sec. 12.11 

In the following, in order to use a shorthand notation, we will call processes such as the 
one depicted in Fig.[T] "EW Vjj production", or VBF production of W/Z plus two jets, since 
we consider these processes with the kinematic cuts typical for the selection of VBF (see 
Sec. E)). It should be understood that, in spite of this notation, multiperipheral diagrams 
like (e) and (f) are included, even though they cannot be represented as the production of a 
weak boson, followed by its decay into two leptons. 

The structure of the paper is as follows: in Sec. |H we outline the calculation of the tree- 
level diagrams, of real-emission contributions and of the virtual corrections. We dedicate 
Sec. 12.31 to the discussion of the virtual contributions, with some of the analytical details 
relegated to Appendix ^ A list of checks which we have performed on our calculation 
concludes Sec. El Additional features of our Monte Carlo program, like the gauge invariant 
handling of finite W and Z widths, the inclusion of anomalous WW'j and WWZ couplings, 
the approximations with regard to crossed diagrams in the presence of identical quark flavors, 
the singularities for incoming photons and the choice of parameters, will be discussed in 
Sec. El We then use this Monte Carlo program to present first results for EW Vjj production 
at the LHC. Of particular concern is the scale dependence of the NLO results, which provides 
an estimate for the residual theoretical error of our cross-section calculations. We discuss 
the scale dependence and the size of the radiative corrections for various distributions in 
Sec. m Conclusions are given in Sec. 
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Figure 1: Feynman graphs contributing to the process uc^ dci^v^ at tree level. For 
the generic VBF process discussed in this paper, seven Feynman- graph topologies 
contribute at tree level: the six topologies shown plus an additional bremsstrahlung 
graph, with the vector boson emitted off the final-state charm quark [mirror image 
of graph (b)j. 



2 Elements of the calculation 



The structure of the three processes under consideration, pp^l^UijjX, pp-^i'UijjX 
and pp^£~^£~jjX, is very similar. A discussion of any single one of them is sufficient 
to clarify our procedures for all, and we use W'^ production, i.e., the calculation of the 
pp — > I'^v^jjX cross section, for this purpose. Mutatis mutandis, all the considerations apply 
to the other processes too. 



4 



2.1 Approximations and general framework 

At tree level, the topological structure of the generic subprocesses contributing to EW Wjj 
production is depicted in Fig. ^ Two additional classes of diagrams appear in case of 
identical quark flavors on two of the fermion lines: 

- diagrams where both the two virtual vector bosons are time-like. They correspond 
to diagrams called conversion, abelian and non-abelian annihilation in Ref. and 
contain vector-boson pair production with subsequent decay of one of the weak bosons 
to a pair of jets. Pars pro toto, we call this class vector-boson pair production in the 
following. 

- diagrams obtained by interchange of identical initial- or final-state (anti)quarks, such 
as in the uu dul^vi or du — > ddl^vi subprocesses. 

These additional diagrams are obtained from the ones shown in Fig. Q by crossing. In our 
calculation, we have neglected contributions from vector-boson pair production completely. 
In addition, any interference effects of the second class with the graphs of Fig. ^ are ne- 
glected. This is justified because, in the phase-space region where VBF can be observed 
experimentally, with widely-separated quark jets of very large invariant mass, the neglected 
terms are strongly suppressed by large momentum transfer in one or more weak-boson pro- 
pagators. Color suppression further reduces any interference terms. We have checked with 
MadEvent ^7] that, at LO, the diagrams that we have not considered and interference effects 
contribute less than 0.3% to our final results in e.g. Fig. HI Since we expect QCD corrections 
to the neglected terms to be modest, the above approximations are fully justified within the 
accuracy of our NLO calculation. 

Fermion masses are set to zero throughout, because observation of either leptons or (light) 
quarks in a hadron-coUider environment requires large transverse momenta and hence sizable 
scattering angles and relativistic energies. For the t-channel processes which we include, we 
have used a diagonal form (equal to the identity matrix) for the Cabibbo-Kobayashi-Maskawa 
matrix, Vckm- This approximation is not a limitation of our calculation. As long as no final- 
state quark flavor is tagged (no c tagging is done, for example), the sum over all flavors, 
using the exact Vckm, is equivalent to our results, due to the unitarity of the Vckm matrix. 

2.2 Tree-level diagrams and real corrections 

For the Wjj Born amplitude, we need to add the contributions from the 10 Feynman graphs 
shown in Fig. n](Z and 7 propagators counted as different diagrams), and sum cross sections 
of all subprocesses producing plus two jets. The same is true for W~ production. For 
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Figure 2: Examples of Feynman amplitudes with an initial gluon. Graphs like (a) 
and (h), with the gluon coupled to the initial quark line, correspond to vector-boson 
pair production and are eliminated. The two gauge-invariant subsets of graphs 
like (c) and (d), with the gluon coupled to the final-state quark pair, contain all 
g^qq splitting contributions and are included in our calculation. 



the case of Zjj production, amplitudes which correspond to neutral-current exchange (no 
change of quark flavors) receive contributions from 24 Feynman graphs at tree level. To 
obtain the real-emission diagrams, with a final-state gluon, one needs to attach the gluon 
to the quark lines in all possible ways. For the diagrams in Fig. Q this gives rise to 45 
real-emission graphs. 112 different Feynman graphs contribute to real-emission corrections 
to Zjj production via neutral-current exchange. 

The contributions with an initial-state gluon are obtained by crossing the previous dia- 
grams, promoting the final-state gluon as incoming parton, and an initial-state (anti)quark 
as final-state particle. We again remove all diagrams where two time-like, final-state vector 
bosons appear such as gu-^ l^p^dZ* , with Z* ^ cc. Such diagrams, for consistency, must 
be removed since we have not considered the corresponding Born contributions. Figure |21 
clarifies this issue: we drop all initial-gluon contributions in which the gluon couples to the 
fermion line of the initial quark or antiquark. In fact, these diagrams are strongly suppressed 
when VBF cuts (see Sec. EJ are applied to the final-state jets. 

Our Monte Carlo program computes all amplitudes numerically, using the formalism of 
Ref. The Born amplitudes for W and Z production are taken from Ref. [^. The 
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real-emission amplitudes for Z production were first given in Ref. (?]. The corresponding 
amplitudes for W production were partially programed at the time. We have finalized and 
tested them for the present application. 



2.3 Virtual corrections 

At NLO, we have to deal with soft and coUinear singularities in the virtual and real-emission 
contributions. Our calculation uses the subtraction method of Catani and Seymour ^H] to 
cancel the soft and collinear divergences between virtual and real-emission diagrams. Since 
these divergences only depend on the color structure of the external partons, the subtraction 
terms encountered for EW V production are identical in form to those found for Higgs 
boson production in VBF. Thus, we can use the results described in Ref. ^H] for the case 
at hand. The main difference is that the finite parts of the virtual corrections are more 
complicated than for Hjj production (where only vertex corrections were present). 

The QCD corrections to EW V production appear as two gauge-invariant subsets, 
corresponding to corrections to the upper and lower fermion lines in Fig. ^ Due to the 
color singlet nature of the exchanged electroweak bosons, any interference terms between 
subamplitudes with gluons attached to both the upper and the lower quark lines vanish 
identically at order a^. Hence, it is sufficient to consider radiative corrections to a single 
quark line only, which we here take as the upper one. Corrections to the lower fermion line 
are an exact copy. 

In computing the virtual corrections, we have used the dimensional reduction scheme [20] '■ 
we have performed the Passarino-Veltman reduction of the tensor integrals in c? = 4 — 2e 
dimensions, while the algebra of the Dirac gamma matrices, of the external momenta and of 
the polarization vectors has been performed in = 4 dimensions. 

We split the virtual corrections into two classes: the virtual corrections along a quark 
line with only one weak boson attached and the virtual corrections along a quark line with 
two weak bosons attached. 

I. The virtual NLO QCD contribution to any tree level Feynman subamplitude m!^^ 
which has a single electroweak boson V (of momentum q) attached to the upper fermion 
line. 



appears in the form of a vertex correction, which is factorisable in terms of the original Born 
subamplitude 



g(fci)^g(fc2) + \^(g). 



(2.1) 



M^) = M 




r(i + e) 



■ 2 



3 



e 



+ Cvirt + O (e) 



(2.2) 
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Figure 3: Virtual corrections for a fermion line with two attached electroweak 
bosons, Vi{qi) and V2{q2)- The finite part of the sum of these graphs defines the 
reduced amplitude Air {qi,q2) of Eq. \2.5]) . 



Here is the renormalization scale, and the boson virtuahty = —{ki — A;2)^ = —q"^ is the 
only relevant scale in the process, since the quarks are assumed to be massless, k\ = = 0. 
In dimensional reduction, the finite contribution is given by Cvirt = — 7 (cvirt = — 8 
in conventional dimensional regularization). 

II. The second class of diagrams are the virtual QCD corrections to the Feynman graphs 
where two electroweak bosons Vi and V2 (of outgoing momenta qi and ^2) are attached to 
the same fermion line (see, for example, the upper quark line in Fig. ^ (a, b)). It suffices 
to consider one of the two possible permutations of Vi and V2, as depicted in Fig. IHl The 
kinematics is given by 

q{k^) q{k2) + ^i(gi) + ^2(^2) , (2.3) 

where kf = k2 = and momentum conservation reads ki = k2 + qi + q2- In the following, it is 
convenient to use the Mandelstam variables for a 2 ^ 2 process which we take as qq — * V1V2. 
We then define 

s = (fci-fca)' = (gi + g2)' , t = [k^-qif = {k2 + q2f , u = (^1-^2)' = (^z + gi)' . (2.4) 

In order to use the same notation as in Eq. ()2.2j) . we define = 2ki ■ k2 = —s. 

The two electroweak bosons are always virtual in our calculation, i.e., the effective polar- 
ization vectors ei(gi) and t2{q2) actually correspond to fermion currents (the charm-quark 
current and the leptonic-decay currents in the Feynman graphs of Fig. ^ (a, b)). Since 
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fermion masses are neglected, current conservation implies transversity of the effective po- 
larization vectors: ei • gi = £2 ■ q'2 = 0. The expressions that we give in Appendix 1X1 exploit 
this relationship. Our numerical code permits to switch on the missing ei ■ qi and 62 ■ q2 
terms, allowing us to test gauge invariance. Due to the trivial color structure of the corre- 
sponding tree-level diagram, the divergent part (soft and collinear singularities) of the sum 
of the four diagrams in Fig. Elis a multiple of this Born subamplitude, just like for the vertex 
corrections, 



boxline 



+ 



r(i + 6) 



3 

e 



Cvirt 



47r 



(2.5) 



Here r denotes the quark chirality and the electroweak couplings g^-^ follow the notation of 
Ref. jlHI, with, e.g., g]/ = Qf, the fermion electric charge in units of |e|, g^-^ = l/{V2sm6w) 
and g^^ = (T^f — Qf sin^ 9w)/{sm9w cos9w), where 9]^ is the weak mixing angle and T3J is 
the third component of the isospin of the (left-handed) fermions. 



A finite contribution of the virtual diagrams, which is proportional to the Born amplitude 
(the Cvirt term) , is pulled out in correspondence with Eq. ()2.2|) . The remaining non- universal 
term, A^t(q'I) 92)5 is also finite and can be expressed in terms of the finite parts of the 
Passarino-Veltman Bq, Cq and Dij functions, which we denote as Bq, Cq and Dij. Analyt- 
ical expressions for these functions, along with the expression for Air{qi,q2), are given in 
Appendix |X1 

An equivalent form for Eq. ()2.5j) has been derived where all the Dij have been reduced 
to Bq, Cq and Dq functions. We have checked numerically that the two expressions agree 
within the numerical precision of the two FORTRAN codes. 



The factorization of the divergent contributions to the virtual subamplitudes, as multiples 
of Ai^B J implies that the overall infrared and collinear divergence multiplies the complete 
Born amplitude (the sum of the Feynman graphs of Fig. . We can summarize this result for 
the virtual corrections to the upper fermion line by writing the complete virtual amplitude 

My as 



+ 



Arc 



R) 



in 




J^B -. 

Arc 



ra + e) 



e 



Mriqi, q2)g^'^'g^'^' + MAq2, qi)g'^'^''g'^'^' + Oie) 




r(i + e) 



e 



(2.6) 



where Aiv is finite. The interference contribution in the cross-section calculation is then 
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given by 



2Re[MvM] 



Ztx 




— ^ he 



^virt 



2 Re 



MvM% 



B 



(2.7) 

This expression replaces the analogous result for the NLO corrections to qq^qqH, Eq. (2.11) 
in Ref . [12] ■ The divergent piece appears as the same multiple of the Born amplitude squared 
as in the qq-^qqH cross section. It cancels explicitly against the phase-space integral of the 
dipole terms (see Ref. ^ and Eq. (2.10) of Ref. [H]) 



(1(e)) = \M 



R) 



B\ 



2n 



r(i + 6) 



2 3^ 
+ - + 9 



4 



e2 e 



-TT 



(2.8) 



which absorbs the real-emission singularities. After this cancellation, all remaining integrals 
are finite and can, hence, be evaluated in d = 4 dimensions. This means that the values of 
M. B and Aiv need to be computed in 4 dimensions only and we use the amplitude techniques 
of Ref. jlHI to obtain them numerically. 



2.4 Checks 

We have verified, both analytically and numerically, the gauge invariance of Eq. ()2.6p : once 
the extra ei ■ qi and €2 ■ q2 terms have been re-inserted in this expression, the individual 
finite subamplitudes Air{qi,qj) vanish upon the replacements ei^qi or €2— >g2- This is a 
strong check of the tensor reduction and manipulation of the virtual contributions depicted 
in Fig. El 

We have taken the Born amplitudes for W and Z production from Ref. jB] and use the 
real-emission amphtudes of Ref. jTj for Z production. In addition, the Zjj results at the 
Born level were successfully checked with COMPHEP code j2I] • For W production, the real- 
emission amplitudes were obtained by modifying the previously tested Zjjj amplitudes [Zj . 
We have generated equivalent amplitudes with MadGraph [T7], checking their consistency 
numerically. 

For the case, we have built two totally-independent codes. This has allowed us to 
check the overall structure of the dipole-formalism terms (common to all the vector-boson 
fusion processes), and to compare tree- level, real-emission and virtual amplitudes. The two 
codes agree within the numerical precision of the two FORTRAN programs for the total 
cross sections and for final-state kinematic distributions. 
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3 The parton-level Monte Carlo 



The cross-section contributions discussed above have been implemented in a parton-level 
Monte Carlo program for i'^i'ejj, ^~^ejj and production at NLO in QCD, which is 

very similar to the program for Hjj production by weak-boson fusion described in Ref . J3] . 
As in our previous work, the tree-level and the finite parts of the virtual amplitudes are 
calculated numerically, using the helicity-amplitude formalism of Ref. JH] • The Monte Carlo 
integration is performed with a modified version of VEGAS |22]- While many aspects of our 
present calculation are completely analogous to those described in Ref. [1^], several new 
problems appear for the vector-boson production processes which require explanation. 

In order to deal with W/Z boson decay 

W/Z{p,, ^h{p,,) + £2{pe,) , (3.1) 

we have to introduce a finite W/Z width, Ty, in the resonant poles of the s-channel weak- 
boson propagators. However, in the presence of non-resonant graphs, like those of Figs, ^e) 
and (f), this introduces changes in a subclass of Feynman graphs only, which leads to a 
violation of electroweak gauge invariance, which is guaranteed for the zero- width amplitudes. 
Such non-gauge-invariant finite-width effects can lead to huge unphysical enhancements at 
very small photon virtuality and should be avoided For the case at hand, transverse- 
momentum cuts on the two final-state tagging jets (see Sec.|3I) largely eliminate the dangerous 
phase-space regions with low- virtuality gauge bosons. Nevertheless, a careful handling of the 
finite-width effects is called for. 

We have accomplished this using two different schemes. 
I. In the overall- factor scheme [21] , one multiplies all the diagrams shown in Fig. ^ and all 
virtual and real-emission contributions as well, by an overall factor 

{pe, + pe^r - ml + zmyTy ^ ^^''^ 

where Ty has been assumed to be constant. This way, close to resonance [(j9^^ +^^2)^ ~ ""^y]; 
where the sum of the diagrams is dominated by the vector-boson propagator, we recover the 
result of the resonance approximation. Away from resonance, and, thus, in a subdominant 
phase-space region, the error that we make, by multiplying all the diagrams by the factor in 
Eq. ()3.2|) . is of the order of Ty/my ~ 2.7%, for both Z and W boson production. 

The advantage of this scheme is that it preserves full SU{2) x U{1) gauge invariance, 
since the gauge-invariant set of zero-width diagrams is multiplied by an overall factor. 

II. In the complex-mass scheme 25], one globally replaces my my — imyTy, also in 
the definition of the weak mixing angle, siia^ 9w = 1 — m^/m'^^. We have implemented a 
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modified complex- mass scheme where we replace rriy —>■ niy — imv^v in the weak-boson pro- 
pagators appearing in Fig.[TJ but we keep a real value for sin^ 9w- With this prescription, the 
electromagnetic Ward identity relating the tree-level triple-gauge-boson vertex, —ieV^^r^, 
and the inverse W propagator, is preserved 122] 

(gi - q2)^,'^ww^ = - i{Dw)~^p{q2) ■ (3.3) 



This relation removes potential problems with small photon propagators, where gauge- 
invariance- violating terms, proportional to Tw/mw, may be enhanced by factors E^/q'^, 
where the hard scale Et is set by typical transverse momenta of the process. The cor- 
responding enhancement for Z-boson propagators is of order E^/{\q'^\ + m|) and, hence, 
small for the energies available at the LHC Also, we note that the imaginary part of 
sin^ 9\Y = 1 — i^w ~^ '^"^VKriy)/(m| — inizTz), in the full complex-mass scheme, is 200 
times smaller than the real part and hence well below the naive expectation Tv/fny ~ 2.7% 
for the size of finite-width corrections. 

We have used the two different schemes to compute total cross sections with VBF cuts 
and find agreement at the level of the 0.5% or better. This ambiguity, thus, represents a 
minor contribution to higher-order electroweak corrections. 

Inspection of the Feynman graphs of Fig.^shows that the non-abelian triple-gauge-boson 
vertices (TGV) enter via the WWZ and WW'-f couplings in diagrams like Fig. ^(d). These 
graphs receive QCD vertex corrections only and, therefore, factorize according to Eq. ()2.2|) 
in terms of the tree-level TGV graphs, independent of the form of the TGV. In particular, 
the presence of anomalous WWZ or WWj couplings can easily be taken into account by a 
simple modification of the Born amplitude. Our program supports anomalous couplings k^, 
Kz, X-y, Xz etc. j2Zj and thus allows to extend the analysis of anomalous- coupling effects in 
vector-boson fusion processes 0j to NLO QCD accuracy. 

The requirement of two observable jets, of finite transverse momentum (see Sec. |3I), is 
sufficient to render the LO cross section for EW Wjj and Zjj events finite. At NLO, initial- 
state coUinear singularities appear. For g^qq and q—>-qg splitting, these are properly taken 
into account via the renormalization of quark and gluon distribution functions. An additional 
coUinear divergence exists, however, because of the presence of t-channel photons in tree- 
level graphs, such as in Fig. ^(a, b, d, e). Real-emission corrections lead to Feynman graphs 
such as the one shown in Fig. |21 (d): the final-state d and u quarks may lead to observable 
jets, allowing vanishing momentum transfer for the virtual photon and a corresponding 
coUinear singularity, representing, in the case shown, a QED correction to the LO process 
g'~f—>-duW^. This singularity would have to be absorbed into the renormalization of the 
photon distribution function inside the proton. Alternatively, one may impose a cut, |t| > 
Q7,min' '^^ ^hc virtuality of the photon and replace the missing piece by the p'j—>VjjX cross 
section, folded with the appropriate photon density in the proton j211l2Hl. We have chosen 
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this latter approach: all divergent amplitudes are set to zero below = 4 GeV^ and 

py—>-VjjX is considered to be a separate electroweak contribution to Vjj events, which we 
do not calculate here. 

When imposing typical VBF cuts, with their large-rapidity separation and concomitant 
invariant mass of the two tagging jets, the p'j^VjjX contribution to the EW Vjj cross 
section is quite small. For the VBF cuts defined in the next section, with pxj > 20 GeV and 
a rapidity separation of the two tagging jets of Ai/jj > 4, the NLO W^jj cross section, for 
example, increases by a mere 0.2% when lowering the photon cutoff to Q^ min = 0.1 GeV^ 
from our 4 GeV^ default value^. This number increases to 0.7% for ^Vjj > 2. Because 
these contributions are negligible, we have not yet implemented the calculation of this small 
missing piece in our program. 

In the computation of cross sections and distributions presented below, we have used 
the CTEQ6M parton distribution functions (PDFs) [SHI with a,(mz) = 0.118 for all NLO 
results and CTEQ6L1 parton distributions for all LO cross sections. The CTEQ6 fits include 
b quarks as an active fiavor. For consistency, the b quark is included as an initial- and/or 
final-state massless parton in all neutral-current processes, i.e., we include only those pro- 
cesses with external b quarks, where no internal top-quark propagator appears via the btW 
vertex, being forbidden by Feynman rules. Top-quark contributions, obviously, go beyond 
our massless-fermion approximation and would have to be treated as a separate process. 
Allowed neutral-current processes with b quarks appear for Z production only. The 6-quark 
contributions are quite small, however, affecting the Z-boson production cross section at the 
1% level only. 

We choose mz = 91.188 GeV, rriw = 80.419 GeV and the measured value of Gp as our 
electroweak input parameters, from which we obtain agED = 1/132.51 and sin^ 9w = 0.2223, 
using LO electroweak relations. The decay widths are then calculated as Tw = 2.099 GeV 
and Tz = 2.510 GeV, which agrees with their Particle Data Group values at the level of 
0.9% and 0.6% respectively, which is better than the overall theoretical uncertainty we are 
striving for. 

In order to reconstruct jets from the final-state partons, the kj- algorithm ^T], as de- 
scribed in Ref. jS2|, is used, with resolution parameter D = 0.8. 

^The finite proton mass provides an absolute lower bound on the photon virtuality, ^ mp{myjj/xs)'^, 
where mvjj is the invariant mass of the produced system and x denotes the Feynman x of the colored parton 
in the subprocesses for pj^VjjX. We have chosen the lower cutoff of ^^^^^ = 0.1 GeV^ for a very rough 
simulation of the resulting finite photon flux. 
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4 Results for the LHC 



The parton-level Monte Carlo program described in the previous section has been used 
to determine the size of the NLO QCD corrections to EW Vjj cross sections at the LHC. 
Using the fcr algorithm, we calculate the partonic cross sections for events with at least two 
hard jets, which are required to have 

Pt, >20GeV, |%|<4.5. (4.1) 

Here i/j denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four- vector sum of massless partons of pseudorapidity |?7| < 5. The two reconstructed jets of 
highest transverse momentum are called "tagging jets" and are identified with the final-state 
quarks which are characteristic for vector-boson fusion processes. 

We consider decays Z^£~^i~ and W^iu^ into a single generation of leptons. In order to 
ensure that the charged leptons are well observable, we impose the lepton cuts 

PTi > 20 GeV , \rii\ < 2.5 , ARji > 0.4 , (4.2) 

where Rje denotes the jet-lepton separation in the rapidity-azimuthal angle plane. In ad- 
dition, the charged leptons are required to fall between the rapidities of the two tagging 
jets, 

yj,min ^ yj,max ■ (4-3) 

We do not specifically require the two tagging jets to reside in opposite detector hemi- 
spheres for the present analysis. Backgrounds to VBF are significantly suppressed by re- 
quiring a large rapidity separation of the two tagging jets. Unless stated otherwise, we 
require 

^Vjj = IVn - Vhl > 4= ■ (4-4) 

Cross sections, within the cuts of Eqs. (j4.H) - (j4.4|) . are shown in Fig. |1[ for Wjj produc- 
tion, and in Fig. El for the Zjj case. In both figures, the scale dependence of the LO and 
NLO cross sections is shown for fixed renormalization and factorization scales, /i/j and /i^, 
which are tied to the masses of the produced vector bosons my 

(J'R = ^Rmv, fXF = ^F'mv. (4.5) 

The LO cross sections only depend on fip = ^my. At NLO we show three cases: (a) 
= in = i (red solid line); (b) = = ^ (blue dot-dashed line); and (c) = ^, 

^i? = 1 (green dashed line). While the factorization-scale dependence of the LO result is 
sizable, the NLO cross sections are quite insensitive to scale variations: allowing a factor 2 
variation in either directions, i.e., considering the range 0.5 < < 2, the NLO cross sections 
change by less than 1% in all cases. 
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Figure 4: Scale dependence of the total cross section at LO and NLO within the cuts 
of Eqs. for W~ and W'^ production at the LHC. The decay branching 

ratio of the W is included in the definition of the cross section, here and in all 
subsequent figures. The factorization scale fip and/or the renormalization scale 
have been taken as multiples of the vector-boson mass, ^mw , and ^ is varied in the 
range 0.1 < ^ < 10. The NLO curves are for ^ip = ^^r = (solid red line), 

Hf = fnw CLf^d hb. = £,mw (dashed green line) and fiR = mw and variable 
(dot-dashed blue line). The dotted black curve shows the dependence of the LO cross 
section on the factorization scale. At this order, as{nR) does not enter. 



As a second option, we have considered scales tied to the virtuahty of the exchanged 
electroweak bosons. Specifically, independent scales Qi are determined as in Eqs. ()2.2|1 
and ()2.5j) for radiative corrections on the upper and on the lower quark line, and we set 

^^Fi = ipQi , l^m = inQi ■ (4.6) 

This choice is motivated by the picture of VBF as two independent deep-inelastic scattering 
type events, with independent radiative corrections on the two electroweak-boson vertices. 
Resulting Vjj cross sections at NLO are about 1% lower for np = = Qi than for 
/^F = A'-R = fT^v- In the following, we refer to the latter choice as the "M scheme" while the 
choice fip = fiR = Qi is called the "Q scheme". As we will see below, a residual NLO scale 
dependence of about l%-2% is also typical for distributions, resulting in very stable NLO 
predictions for Vjj cross sections. 

In addition to these quite small scale uncertainties, we have estimated the error of the 
W^jj cross sections due to uncertainties in the determination of the PDFs. This error is 
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Figure 5: Same as Fig. ^ but for Z production at the LHC, with the Z^fi^fi~ 
branching ratio included in the definition of the cross section, here and in all subse- 
quent figures. 



determined by calculating the total Wjj cross section, within the cuts of Eqs. (j4.ip - (|4.4j) . 
using two different sets of PDFs with errors, computed by the CTEQ [^H] and MRST '3'^ 
Collaborations. Together with the PDF that gives the best fit to the data, the CTEQ6M set 
provides 40 PDFs, and the MRST2001E 30 PDFs, which correspond to extremal plus- minus 
variations in the directions of the error eigenvectors of the Hessian, in the space of the fitting 
parameters. To be on the conservative side, we have added the maximum deviations for each 
error eigenvector in quadrature, and we have found a total PDF uncertainty of ±4% with 
the CTEQ PDFs, and of roughly ±2% with the MRST set. 

For precise comparisons with future LHC data, the residual theoretical error on the jet 
and lepton distributions must be estimated. As a first example, we show the transverse- 
momentum distribution of the highest-p^ tagging jet for W^jj production in Fig. El (a): the 
shape of the pt distribution is fairly similar at LO (red dashed curve) and NLO (black solid 
line). Both curves were obtained with a scale choice of hr = fip = mw- In the right-hand 
panel their ratio to the NLO curve with fin = fip = Qi is shown. The ratio of the two NLO 
distributions deviates from unity by 2% or less over the entire range, which, again, points to 
the small QCD dependence of our calculation. 

In contrast to the stability of the NLO result, the LO curves depend appreciably on the 
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Figure 6: Transverse-momentum distribution of the highest-px tagging jet in 
production at the LHC. In panel (a) the NLO result (solid black line) and the LO 
curve (dashed red line) are shown for the scale choice = l^R = (M scheme). 
In panel (b), we show the ratios of the NLO differential cross section in the M 
scheme (solid black line), of the LO one in the M scheme (dashed red line) and of 
the LO one in the Q scheme (blue dotted line) to the NLO distribution in the Q 
scheme, which is defined via the scale choice fip = fJ-R = Qi- 

scale choice. The blue dotted line and the red dashed line in Fig. IHl (b) give the ratio of 
the LO curves for fip = Qi and fip = mw, respectively, to the NLO result. The shape of 
the LO curves, in particular for a constant scale choice like fip = is quite different 

from the more reliable NLO result. For transverse-momentum distributions we generally 
find that the "dynamical" scale choice = Qi, at LO, better reproduces the shape of 
the NLO distributions, and is thus preferable to a fixed scale. At NLO, or higher order, 
where the definition of the momentum transfer Qi becomes more problematic, the fixed-scale 
choice becomes more natural. However, because of the greater stability of the cross-section 
prediction, the scale selection also becomes less of a phenomenological issue. 

Rapidity distributions of the two tagging jets are shown in Fig. [71 at LO and NLO, and 
for two choices of the rapidity-gap requirement, Ayjj > 2 and Ayjj > 4. The shapes of 
the rapidity distributions for the more central tagging jet, panel (a), and the more forward 
tagging jet, panel (b), are quite similar at LO and NLO. In fact, the K factors for these 
distributions are fairly fiat, and adequately described by a constant value of about 1.1. 
The results in Fig. [7| were obtained for a fixed scale = jJ^R = f^w and are for W~jj 
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Figure 7: W~ production cross section as a function of (a) the smaller and (h) 
the larger absolute value of the two tagging-jet rapidities. Results are shown for a 
rapidity separation between the two tagging jets greater than 2 and 4 (higher and 
lower pairs of curves, respectively). The LO cross section is always slightly below the 
NLO result. Due to the rapidity cut of Eq. J^. the distributions are truncated at 
|y,|=4.5. 



production. Curves for the W^jj and Zjj cross sections are very similar in shape and show 
the preservation of shape between LO and NLO curves. 

While tagging-jet distributions are quite similar for electroweak Wjj and Zjj events at 
the LHC, the presence of two charged leptons in the Zjj case results in somewhat more 
noticeable differences. When considering changes in the lepton cut of Eq. ()4.2|) . the 
transverse momentum of the softer lepton is critical for Z production, while the single charged 
lepton must be considered for Wjj events. These distributions are shown in Fig. |H1 for 
production (top panels) and Z production (bottom panels). At NLO the scale variations are 
again very small, at the 1% level, as demonstrated by the ratios of the NLO pt distributions 
for fJ,F = fi'R = fny and fJ^F = fJ^R = Qi (solid black lines) in Fig. |H1 (b, d). Varying either 
scales by a factor of 2 leads to the same conclusion of l%-2% scale uncertainties for the 
NLO results. Comparing the LO predictions (dashed and dot-dashed curves) with the very 
precise NLO results shows theoretical errors of the order of 10%. Again, as for the jet pr 
distributions discussed earher, the choice fip = Qi is better for simulating the shape of the 
lepton Pt distribution at LO. A fixed scale, fip = f^v, predicts too steep a fall-off at large 
Pt. One should note, however, that for the electroweak Vjj processes considered here, these 
differences are exceptionally small already at LO: the differences between the LO curves in 
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Figure 8: Transverse-momentum distributions of the charged final-state lepton in 
production [panels (a) and (h)] and of the softest of the two final- state leptons in 
Z production [panels (c) and (d)]. The solid black curves in panel (a) and (c) rep- 
resent the NLO cross sections and the red dashed curves the LO ones, for scales 
fJ'R = fJ-F = T^v scheme). Panels (b) and (d) show the ratio of the NLO 
transverse-momentum distribution computed in the M and Q scheme (black solid 
line), and the K factors in the Q (red dashed line) and M (blue dot-dashed line) 
schemes. 
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Fig. ISlare of the order of 10% only. 



In contrast to the lepton transverse-momentum distributions described above, the shape 
of the lepton-rapidity distributions is virtually unaffected by the NLO corrections: an overall 
constant K factor is sufficient to describe NLO effects. Larger changes are found when 
considering angular correlations of the leptons and jets, which we show for Zjj production 
in Fig. ini The top panels show the minimal rapidity between any of the two leptons and 
the two tagging jets, AyJ^™;. As before, the tagging jets are taken as the two highest 
transverse- momentum jets in the event {px selection). The two bottom panels show the 
minimal separation in the rapidity-azimuthal angle plane of the two leptons from any jet 
(not necessarily the two tagging jets) in the event, In both cases, the two scale choices 

for the NLO result show excellent agreement (black solid lines in Fig. El (b, d)). However, 
the dynamical K factors 

K{x) = ^yo/dx 
daLo/dx 

for X = Ay'^'^i and x = -R™;™ show qualitatively different behavior. While K{Ay^^i) is fairly 
constant, i.e., the shape of the distribution is well described by the LO approximation, the 
minimal lepton-jet separation, da/dR^l'^, shifts noticeably to smaller values at NLO. This 
behavior was to be expected, since additional parton emission in the higher-order calculation 
reduces lepton isolation. What is remarkable, then, is that the selection of the tagging jets 
as the two highest-pr jets does not affect the lepton-tagging jet separation. As for the Higgs 
boson case JHI; this selection of the tagging jets provides excellent correspondence of the 
LO- and NLO-event topology. 

In order to stress this point we show dijet invariant-mass distributions for the recon- 
structed jets (not necessarily the two tagging jets) for W^jj events at LO (red dashed lines) 
and at NLO (solid black lines) in Fig. ^1 The distribution with respect to the minimal dijet 
invariant mass in the event is shown in Fig. (a) while Fig. (b) uses the invariant mass 
of the two tagging jets, mtags- At LO, there are only two final-state quarks of pt > 20 GeV 
in each event and, hence, the two curves are identical. At NLO, additional parton emission 
provides for soft third jets which form low invariant- mass pairs with one of the tagging jets, 
and this pair shows up as a low-mass peak in dcr/fim™™. Generic selections of the two tag- 
ging jets in a multijet environment tend to pick up some of these low-mass pairs and lead 
to substantial differences in the invariant-mass distribution of the two tagging jets at LO 
and at NLO. The pt selection of tagging jets, which we have used throughout and for which 
results are shown in Fig. (b), is remarkable in that it preserves the shape of the tagging 
jet invariant-mass distribution, dcr/iimtags, when going from LO to NLO. 
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Figure 9: Angular correlations of leptons and jets in Z production. Panels (a) 
and (h) show the minimum rapidity separation between the two leptons and the two 
tagging jets. Panels (c) and (d) are for the minimum rapidity- azimuthal angle sepa- 
rations between the leptons and any reconstructed jets (not necessarily the two tagging 
jets). The NLO differential cross sections are shown in black solid lines, while the 
LO ones are displayed as red dashed lines. Scales are fixed in the M scheme. Pan- 
els (h) and (d) show the ratio between the two NLO differential cross sections in the 
M and Q scheme (solid black lines) and their respective K factors. 
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Figure 10: Dijet invariant-mass distributions for production, with scales in 
the M scheme. Shown are (a) the minimum dijet invariant-mass distribution for 
any final-state reconstructed jets (not necessarily the two tagging jets) and (b) the 
invariant mass of the two tagging jets. NLO results are shown in solid black lines, 
while the red dashed lines are for LO distributions. 

5 Conclusions 

Vector-boson fusion at the LHC represents a class of electroweak processes which are 
under excellent control perturbatively. This has been known for some time for the most 
interesting process in this class: Higgs boson production via VBF has a modest K factor 
of about 1.05 for the inclusive production cross section ^3] and this result also holds when 
applying realistic acceptance cuts 

In the present paper, we have extended this result to the electroweak production of 
W and Z plus two jets, when the final-state particles are in a kinematic configuration 
typical of VBF events. More precisely, we have calculated the NLO QCD corrections to 
electroweak production of iiyejj and at LHC, and we have implemented them in a 

fully-flexible NLO Monte Carlo program. K factors are of the same size as for the Higgs 
boson production process, typically ranging between 1.0 and 1.1 for most distributions. What 
is more important is the stability of the NLO result: residual scale dependence is at the 2% 
level or below. This is smaller than the present parton-distribution-function uncertainties, 
which we have calculated for the W^jj cross sections. We estimate 4% PDF errors using 
CTEQ6M parton distributions and roughly half that size using MRST2001E PDFs. 
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Given the excellent theoretical control which we now have for EW Vjj production, these 
processes can be used as testing grounds for Higgs boson production in VBF: techniques 
should be developed to measure hadronic properties, like forward-jet tagging efficiencies or 
central-jet- veto probabilities, in Wjj or Zjj production at the LHC and to extrapolate 
these results to Higgs boson production, thus reducing the systematic errors for Higgs boson 
coupling measurements. We leave such applications for the future. 



Acknowledgments 

Part of this work was done at LAPTH in Annecy and D.Z. would like to thank the 
members of the laboratoire for their hospitality. This research was supported in part by the 
University of Wisconsin Research Committee with funds granted by the Wisconsin Alumni 
Research Foundation and in part by the U.S. Department of Energy under Contract No. DE- 
FG02-95ER40896. CO. thanks the UK Particle Physics and Astronomy Research Council 
for supporting his research. 



A Virtual corrections 



In this appendix, we give the expression for the finite, reduced amplitude Air{qi, Q'2) that 
appears in Eqs. (j2.5p and (j2.6p . in terms of Bq, Cq and Dij functions. Here Bq, Cq and Dij 
are the finite parts of the Passarino-Veltman Bq, Cq and Dij functions and are given 
explicitly below. We have also derived A^t-(5'i,?2) in terms of Bq, Cq and Dq functions, but 
do not show this expression here, due to its length. We write 

Mriqi, q2) = Hk2) [ci^, + C2^ + (j^i - Q + Cb^2 ih + Q A] ^-^^ih) , (A.l) 

where ei = ei(gi) and €2 = ^2(^2) are the effective polarization vectors of the two electroweak 
gauge bosons. The coefficient function ci = ci(gi, ^2) is given by 

Ci = 2e2- k2T,{q2,tj -2 Di2{k2,q2,qi) + D2A{k2,q2,qi) ^2 ■ k2 {^qf + q^ - Ss - ' 

- 2 Di2{k2,q2,qi) - D2i{k2,q2,qi) £2 ■ gi (^2 - ^) 

+ 4 -Du{k2, q2, qi)e2 ■ k2S - -012(^2, ^2, gi)e2 ■ kit + Di3{k2, q2, qi)(^2 ■ k2 (qj - s - t 

+ Di3{k2, q2, qi)t2 ■ qiql - -D2i(/^2, ^2, qi)e2 ■ k2S - D22{k2, q2, gi)e2 ■ k2t 

- -^22(^2, g2, gi)e2 ■ qiql + D22.{k2, q2, qi)e2 ■ k2qf + D2b{k2, q2, gi)e2 ■ k2 (ql - s - 2t 

- D2e{k2, g2, gi)e2 ■ /c2 (^2 - « - ^) + D26{k2, g2, gi)e2 ■ qit + 2D27{k2, g2, qi)^2 ■ qi 
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- D^2{k2, q2, qi)e2 ■ hql - D^i{k2, q2, qi)e2 ■ k2{ql - t) 
+ D36{k2, q2, qi)e2 ■ k2 (2g| -tj + Dzj{k2, q2, qi)e2 ■ k2qf 

+ -035(^2, ?2, qi)e2 ■k2(ql - s -t) + D^^{k2, ?2, qi)^2 ■ k2 [ql + ql - 

- Dsg{k2, q2, qi)e2 ■ k2ql - D3io{k2, ^2, 9i)e2 ■ k2 (q^ + 2ql - 2s - tj 

- 4^3ii(/c2, ?2, qi)e2 ■ k2 + GDsi2{k2, q2, qi)e2 • k2 + 2Dzi^{k2, ?2, qi)^2 ■ qi 



where 



Te(q',t) 



t-q' 



2t + 
t-q^ 



+ 2Bo{q') + l-2q'Co{q\t) 



is defined in terms of the finite parts of the Bq and Cq functions 



So(g2) = 2 - In 



q^ + iO+ 



and 



Coiq',t) 



1 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



2{t - V s s 

These expressions are obtained by pulhng a common factor r(l + e)(— s)~'^ = r(l + e)/{Q^y 
out of all amplitudes and Passarino-Veltman functions, e.g., 

d^k 1 r(l + e)r(l-e)2 



Bo{q') 



r(i + e 



e r(2 - 2e) 



io+)- 



l + 2-ln^ + C.(.) 
e s 



r(i + 6) 



1 



+ Bo{q^) + O{t) 



(A.6) 



For the other coefficient functions q = Ci(gi, ^2) we find 

C2 = -2 -Di2(A;2,?2,gi) + -D24(A;2,?2,gi) ei- k2{ql + ql- s -2t}j + ei- q2{ql- s -?,t^ 

+ 4 Di3{k2, q2, gi)ei ■ A;2gi - ^i3(/^2, 92, gi)ei ■ /i:i(2s + t) + £"22(^2, ^2, gi)ei ■ fci'?2 

- -023(^2, g2, gi)ei • k2t + D2z{k2, q2, qi)ei ■ q2 (qj - t) - D24{k2, q2, qi)ei ■ kiql 

+ D25{k2, 92, gi)ei ■ k2qj + 1)25(^2, ^2, gi)ei ■ ki (^qj - 2s - t) + 1)26(^2, ^2, ■ ^2^ 

- D2e{k2, q2, qi)ei ■ ki [ql - sj - 2D27{k2, q2, qi)ei ■ q2 + D^2,{k2, q2, qi)ei ■ k2ql 

+ D3s{k2, 92, 9i)ei • q2qj + -037(^2, 92, 9i)ei • ki (^qj - s - tj + D3sik2, 92, 9i)ei ■ /i;i92 

- -^39(^2, 92, 9i)ei • /ci (91 + 92 - - -^310(^2, 92, 9i)ei • /^i (92 - ^) 

+ 2L>3ii(/c2, 92, 9i)ei • /C2 + 2D3i2(A;2, 92, 9i)ei • 92 - 6-D3i3(A;2, 92, 9i)ei • ki 

+ 2e,-kiT,(ql,t) , (A.7) 
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Cq = [-Di2(A;2, g2, Qi) + -^24(^2, g2, gi)] ei • + 2[4L>i2(A;2, ^2, gi)e2 • ^2^1 • A;2 

+ 3L'i2(A;2, ^2, gi)e2 • ^261 • q2 + -Di2(A;2, ^2, gi)e2 • giei • A;2 - 4/)i3(A;2, ^2, gi)e2 • k2ei ■ k2 

- 2Dis{k2, ?2, ?i)e2 • A;2ei • q2 - 2bi^{k2, q2, qi)^2 ■ qi^i ■ k2 - Di^{k2, ?2, • £2* 

+ 2£»22(fc2, 52, gi)e2 • ^2^1 • g2 - -022(^2, ^2, gO^l • 62^ - 2£»23(/e2, ?2, ?l)e2 ' ^iCl ' ^2 

- 2D23{k2, q2, gi)e2 ■ gi^i • ^2 - D23{k2, q2, qi)ei ■ e2t + 6D24{k2, q2, qi)e2 ■ k2ei ■ k2 

+ 3^24(^2, q2, qi)^2 ■ k2ei ■ q2 + D2i{k2, ^2, qi)^2 ■ qi^i ■ k2 - 6-D25(/c2, ^2, qi)^2 ■ k2ei ■ k2 

- 2L'25(A;2, 92, qi)e2 • ^^261 • q2 - 2L)25(^2, 92, gi)e2 • gi^i • A;2 - D25ik2, q2, qi)ei • €28 

- 4^26(^2, q2, qi)^2 ■ k2ei ■ q2 + 4:D26{k2, q2, qi)e2 ■ qiei ■ k2 + 2ID2&{k2, q2, qi)t2 ■ qi^i ■ q2 

+ D26{k2, q2, qi)ei • 62 (S + 2t) - D32{k2, q2, qi)ei ■ 62^2 + -033(^2, ^2, qi)^! ' 62^1 

+ 2Du{k2, q2, qi)e2 ■ k2€i ■ k2 - 2^)35 (A;2, ?2, ?i)e2 • k2ei ■ k2 

+ D3Q{k2, q2, qi)ei ■ 62 (gf - t) - 2-037(^2, ^2, qi)e2 ■ qi^i ■ k2 

+ 2r»36(/c2, ?2, qi)e2 ■ k2ei ■ q2 + D3r(k2, ?2, • £2 (gl - « - ^) 

+ 2D3s{k2, 92, gi)e2 • gi^i • ^2 + D38{k2, ^2, gi)ei ■ 62 (g? + 2g| - s) 

- 2D3g{k2, q2, qi)e2 ■ qiei ■ q2 - D3g{k2, q2, qi)ei ■ 62 (2ql + ql - sj 

- 2D3io{k2, q2, qi)e2 ■ ^2^1 • q2 + 2£)3io(A;2, ^2, qi)e2 ■ qi^i ■ /c2 

- D3io{k2, q2, qi)ei ■ 62 (2ql - s - 2t) + 4D3i2{k2, 52, qi)ei ■ £2 

- 4L>3i3(A;2, q2, qi)ei ■ 62] , (A.8) 
Cb = -2{ [-C'36(A;2, q2, qi) + D3r{k2, q2, qi) - 2D3io{k2, q2, qi)] {ql - 

+ D3s{k2, ?2, qi) {qf + 2^2) - -039(^:2, ^2, qi) (2^1 + ^2) } - 2[-C>o(A;2, ^2, qi) 

+ Dii{k2, q2, qi) + Di2{k2, q2, qi) - 2Di3{k2, q2, qi) + D24{k2, q2, qi) - -025(^2, ^2, ^i) 

+ D2%{k2, q2, qi) - -037(^:2, q2, qi) - D3s{k2, ?2, q\) + -039(^:2, ?2, qi) + -C>3io(A;2, ?2, ?i) 

+ 2| [-D22(A;2, 52, 5l) + D23{k2, ^2, ^l) " 2D26{k2, ^2, 5l) ^ " 2Z)27(fe, ^2, 5l) 



+ 



D32(k2, 52, 5l)52 - D33(k2, 52, 5l)5l - 6 (-D3i2(/C2, 52, 5l) - D3i3(k2, 52, 5i)) } 



rfe(g2,i) + r6(g2,i) + So(i) - 5 + 



TT 



(A.9) 



with 



Tb{q',t) = {2q' [Bo{t) - Mq')] + tBo{t) - q^MQ^)} ' '^<fCo{q\t) . (A.IO) 
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For the crossed function J\4{q2,qi), the same expressions as above apply, with the obvious 
interchange qi ^ q2, ei ^ £2, and t^u. 



The finite part of the Dq function is defined by 
1 



Do{k2,q2,qi) 



(A.ll) 



This expression is well defined when all invariants, qf, and t, are space-like. In our 
application, we always have one space-like and one time-like weak boson, i.e., exactly one of 
the two quotients t/qf is positive. In the other quotient simply replace the time-like invariant 
hj t^t + iO+ or qf qf + iO+, as in Eqs. (IA.4j) and (IA.5I1 . 



The remaining finite Dij functions are obtained from the above expressions for the Bo, 
Cq, and Dq functions with the usual Passarino-Veltman recursion relations given in Ref. 
adapted to the Bjorken-Drell metric, qf > for a time-like momentum q^. In these recursion 
relations we need the additional finite Bq and Cq functions 

^o(O) = 0, (A.12) 

C'o(/c2,gi + g2) = ^0(5,0,0) = ^^, (A.13) 

s 

while 

Co{qi,q2) = Co{qlqls) (A.14) 

is the infrared- and ultraviolet-finite Co function for massless internal propagators but with 
nonzero invariants qf, g| and s. 
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